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We propose an extension of the law of corresponding states 
that can be applied to systems - such as colloidal suspensions 
- that have widely different ranges of attractive interactions. 
We argue that, for such systems, the "reduced" second virial 
coefficient is a convenient parameter to quantify the effective 
range of attraction. This procedure allows us to give a simple 
definition of the effective range of attraction of potentials with 
different functional forms. The advantage of the present ap- 
proach is that it allows us to estimate the relative location of 
the liquid-vapor and solid-fluid coexistence curves exclusively 
on basis of the knowledge of the pair-potential. 



Van der Waals's Law of Corresponding States ex- 
presses the fact that there are basic similarities in the 
thermodynamic properties of all simple gases. Its essen- 
tial feature is that if we scale the thermodynamic vari- 
ables that describe an equation of state (temperature, 
pressure and volume) with respect to their values at the 
critical point, all simple fluids obey the same reduced 
equation of state. Pitzer [JIJ has given a molecular inter- 
pretation of the Law of Corresponding States for classical 
monoatomic systems using statistical mechanical argu- 
ments. This proof is restricted to systems for which the 
total intermolecular potential can be written as a sum 
over pair potentials in the form: 



U = J2 € V ( r i,i/ (J ) 



(1) 



The essential assumptions are pairwise additivity and the 
fact that the pair potential can be written as an energy 
parameter e times a function of the reduced distance r/a. 
The Law of Corresponding States follows when we as- 
sume that the pair potentials of all substances to which 
the law applies are conformal. Interaction potentials are 
said to be conformal if their plots can be made to super- 
impose by adjusting the values of e and a. With these 
assumptions, the partition function is of the form: 



Q(N,V,T) 



A 3 



(2) 



where g is the same function for all molecules and de- 
pends only on T* — ksT/e, the reduced temperature, 
and p* — N/Vcr 3 the reduced density |§. It then follows 
that many other thermodynamic properties - in particu- 
lar the pressure - are functions of T* and p* only. 

Unfortunately the interactions between real molecules 
are never truly pairwise additive, nor are the pair po- 
tentials of different molecules conformal. Even for in- 
ert gases the conformality of pair-potentials is only fair. 



Moreover, the importance of three-body forces restricts 
the validity of the assumption of pairwise additivity. 
While only a small family of substances can be de- 
scribed by the original form of the Law of Corresponding 
States, many fluids conform quite accurately to extended 
equations-of-state that involve a third parameter. Thus, 
the compressibility factor z can be expressed as 



z = (3P/p = f(T*,p*,x) 



(3) 



where a; is a third parameter that is usually related to 
some characteristic feature of the phase diagram of a sub- 
stance. At first H the critical compression factor z c was 
used for x, but z c is hard to determine with high ac- 
curacy, and a better choice was sought. The slope of 
the vapor pressure curve (at a reduced temperature of 
T r = 0.7), w, turned out to be a more convenient choice 
for x [fUg]. Various equations of the form z(T* , p* ,uj) 
have been presented |6| , which agree well with the ther- 
modynamic properties for several classes of molecular flu- 
ids. 

In this Communication we focus on the effect of chang- 
ing the range of attractive forces in suspensions of spher- 
ical colloids. As the range of attraction varies indepen- 
dently of the hard-core radius a, the effective interactions 
are clearly not conformal. It is known that the phase be- 
havior of colloidal suspensions depends strongly on the 
range of the attractive interactions. However, at present, 
there is -to our knowledge- no extended law of corre- 
sponding states that allows us to make predictions about 
the phase behavior on the basis of the effective pair po- 
tential alone. In fact, a wide variety of non-conformal 
pair potentials have been used to describe the interac- 
tions between colloids with short ranged attraction. It 
is our aim to formulate an extended law of correspond- 
ing states that allows us to compare such different pair 
potentials. In particular, we have considered the square- 
well model Q], attractive Yukawa potentials 2n-n 
Lennard- Jones type potentials [fiof , the a-Lennard- Jones 
potential used in the description of protein-protein inter- 
actions fll],[l2| , an effective potential reproducing the de- 
pletion attractive forces in colloid-polymer mixtures [fl3| , 
and more complex potentials, which include a repulsive 
barrier, i.e. the effective two-body potential for mix- 
tures of nonadditive asymmetric hard spheres Jl4| ). At 
this stage we limit our analysis to the phase behavior 
around the critical point, but our findings should be gen- 
eralized to densities away from the critical region, in the 
spirit of the extended law of corresponding states. 
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We proceed to calculate the scaling parameters (e, a 
and x), which stem from the knowledge of the inter- 
particle potential alone, without any need for further 
experimental measurement. An obvious choice for the 
length scale a is the effective hard core diameter. Some 
care has to be taken in the calculation of a EFF for 
continuous potentials (such as the Lennard- Jones 2n- 
n). According to the Weeks-Chandler- Andersen (WCA) 
method, we separate the potential into attractive v a tt and 
repulsive v rep parts ]l5| , and calculate the "equivalent" 
hard-core diameter for the repulsive part of the potential 
using the expression suggested by Barker u& : 



KFF 



dr 



1 



(4) 



Two parameters are needed to properly describe the role 
of attractions: an energy scale and a second quantity re- 
lated to the range of attraction. At low temperatures, 
the potential energy per particle in the crystalline phase 
is given by the value of the pair-potential at the nearest- 
neighbor separation, multiplied by the number of neigh- 
bors (and divided by two, to correct for double count- 
ing). This is independent of the functional form of the 
potential. This makes v(r m i n ), the depth of the poten- 
tial well, our natural choice for the energy scale e. The 
third parameter that we use is the reduced second virial 
coefficient, i.e. the second virial coefficient Bi divided by 
the second virial coefficient of hard spheres with a diam- 
eter a EFF . The second virial coefficient B-i can be easily 
calculated once the functional form of the potential has 
been specified: 



Bo = 2tt 



dr 



1 



-v(r)/k B T 



(5) 



and the reduced second virial coefficient £?| is defined as 



b; = 



2ir(a EFF ) 3 /3 



(6) 



Note that all three parameters (<r EFF ,e and £?|) can 
be computed only based on the intermolecular poten- 
tial v(r). In this sense, our approach differs from those 
extended corresponding states laws that use experimen- 
tal data to define appropriate scaling parameters. This 
is particularly useful for the description of colloidal sys- 
tems where the topology of the phase diagram changes as 
the range of the attraction is decreased. For instance, it 
would not be feasible to use the properties at the critical 
point as scaling parameters, as the critical point may be 
experimentally inaccessible for sufficiently short-ranged 
attractions. Our working hypothesis is that for a wide 
range of colloidal materials, the compressibility factor 
z is a function of only three parameters, viz. the re- 
duced temperature T* = ksT /t, the reduced density 
p* = N/V(a EFF ) 3 and the reduced second-virial coef- 
ficient B% 



In the literature onhard particles with short-ranged at- 
traction, it is conventional to express the reduced second- 
virial coefficient in terms of a parameter r that is defined 
through the the following equation |17[: 



B* = 1 - — 
2 At 



(8) 



t is a measure for the temperature - low (high) T cor- 
responds to low (high) t. However, r is not a linear 
function of T. 
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(7) 
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FIG. 1. Stickyness parameter r plotted versus the reduced 
temperature T* — ksT/e for different potentials. Circles: at- 
tractive Yukawa, attractive Yukawa, Squares: Lennard Jones 
2n-2, Triangles: a-Lennard Jones 

In Figure ([!]) we have plotted the stickyness parameter 
r as a function of the reduced temperature for some of the 
cases listed in Table |. In the temperature range stud- 
ied, the stickyness parameter increases almost linearly 
with the reduced temperature [T§| ] . The figure shows an- 
other important feature: if the r — T* curves for two 
different potentials are close at any particular tempera- 
ture, they tend to be close for all temperatures studied. 
Such behavior is an indication that the present scheme 
to compare non-conformal potentials is reasonable. As 
can be seen from the figure (and from Table |), the value 
of r - and therefore that of the reduced second virial 
coefficient - at the critical point is remarkably constant 
(around r w 0.1). This fact had been noted earlier by 
Vliegenthart and Lekkerkerker p9[ . In fact, r hardly 
varies between the limit of extremely narrow attractive 
wells (Baxter's adhesive hard-sphere model [0) and the 
(van der Waals) limit of infinitely long-ranged attractive 
wells. Also in models that are a mixture of the Baxter 
and van der Waals model the value of t at the critical 
point varies only slightly p0[ . 

We mentioned above that the reduced second virial 
coefficient is a measure for the range of the attractive part 
of the potential. To make this statement more precise, 
we have to specify what we mean by the "range" of a 
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potential. Here, we take the following route: there is one 
system for which the range of the attractive potential 
is defined unambiguously, namely hard spheres with a 
square- well attraction 

{oo r < a 
-e a < r < Act . (9) 
Xa <r 

A logical choice for a dimcnsionlcss measure for the 
range of the attractive part of the potential is R = A — 1. 
In the spirit of our extended corresponding states ap- 
proach, we now define the range of an arbitrary attrac- 
tive potential to be equal to the range of that square- well 
potential that yields the same reduced second virial co- 
efficient at the same reduced temperature. The reduced 
second virial coefficient of a square- well potential is given 
by 

B* = l+[X s -l](l-e^ T ') , (10) 
and hence 



4 [A 3 - 1] (e 1 ^* - 1) v ' 

Using this mapping onto the square- well system, we have 
computed the effective range of the attractive part of the 
potential for a number of different potential functions 
that have been used to describe colloidal suspensions or 
globular protein solutions. In general, the range of at- 
traction is still temperature dependent. In Table Q we 
have collected the values of R at the temperature corre- 
sponding to the liquid- vapor critical point. In the same 
table, we also give the value for the "stickyness" param- 
eter r at the critical temperature. 

C 
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FIG. 2. The reduced temperature at the liquid-gas critical 
point T* = ksTc/e plotted versus the range of attractions 
R — A — 1 of the equivalent square well system. As the range 
becomes shorter than the threshold value ~ 0.14, the liq- 
uid-gas transition becomes metastable. Filled Circles: attrac- 
tive Yukawa, Filled Squares: Lennard Jones 2n-n, Filled Tri- 
angles: a-Lennard Jones, Open Squares: Square Well, Open 
Circles: effective Colloid-Colloid interaction 

In Figure (||), we show the relation between T*, the 
reduced critical temperature and R, the range of the at- 
tractive potential. In the temperature range studied, the 
relation between R and T* is surprisingly linear - al- 
though, again, we know that this linear relation cannot 
hold for values of R very close to zero - and obeys the 
simple relation: 

T* m 0.26 + 2. li? . (12) 

The range of the attractive part of the potential deter- 
mines whether a given system can exhibit a stable liquid- 
vapor transition or whether this transition is preempted 
by freezing. The disappearance of the liquid- vapor tran- 
sition in systems with short-ranged attraction was first 
noted in theoretical work by Gast, Russel and Hall 
This work has subsequently been placed on a firmer the- 
oretical footing by Lekkerkerker et al (22). Evidence for 
the disappearing of the liquid-vapor critical point comes 
from both simulation |ll|Jjf and experiment |23]| . All au- 
thors agree that the liquid-vapor transition disappears 
for sufficiently short-ranged attraction. However, esti- 
mates differ for the value of R where this change in the 
phase diagram takes place. Estimates for R vary from 
0.1 to almost 0.4. Part of the reason why the different 
estimates for the critical value of R appear inconsistent is 
that the various authors have studied systems with non- 
conformal interaction potentials and, more importantly, 
have used different definitions for the range. The advan- 
tage of the present approach is that we have a unique way 
to define the range of the attractive potential for widely 
different interaction potentials. When we consider the 
available data for the 2n-n Lennard- Jones potentials, the 
a-Lennard- Jones potential, and the attractive Yukawa 
system, we find that in every case the boundary between 
stable and meta-stable liquid-vapor transitions is located 
within a narrow band between R = 0.13 and R = 0.15. 
In Figure (||) critical points plotted to the right of the 
vertical R « 0.14 band refer to a stable transition, while 
points to the left are metastable. To date no simulation 
has computed the threshold value for square-well parti- 
cles. A rough estimate of R rj 0.25 has been calculated 
using a simple van der Waals model for both the fluid and 
the solid phase [^4| and from a simple cell model with 
some phenomenological character p5[ . For R = 0.85, 
theoretical estimates suggest that in this case the critical 
point is stable [ p6[ . 

The predictive power of our approach based on pair 
potentials is expected to break down when three-body 
interactions become important. We have also tested our 
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theory for more complex pair potentials, which included 
a repulsive barrier [p.4[ . Here too we have found devi- 
ations from extended corresponding states behavior: in 
several cases the calulated r parameter, at the critical 
point, lies much below the constant value of 0.10, and the 
mapping onto the equivalent square-well system yields 
unphysically small attraction ranges R. The repulsive 
barrier alters the effective size of the particle, especially 
at low temperatures, but our WCA decomposition of the 
effective potential limits the repulsive contribution to the 
hard body term. 

In summary, we have formulated a simple extended 
corresponding states principle that allows us to make 
predictions about the topology of the phase diagram of 
suspensions of spherical colloids with variable range at- 
traction. The scaling parameter a, e and 13% can all 
be derived directly from knowledge of the pair-potential. 
Moreover, this procedure allows us to give an unambigu- 
ous definition for the range of the attractive part of the 
potential. By analysing a number of simulation data for 
different model systems, we find that the liquid-vapor 
transition becomes metastable with respect to the freez- 
ing transition when the range of the attraction becomes 
less than approximately 0.14. 
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± c 


T 


R 


Square Well § 


2.61 


0.0765 


1.000 




1.79 


0.0766 


0.750 




1.27 


0.0942 


0.500 




1.01 


0.0924 


0.375 




0.78 


0.1007 


0.250 


Yukawa [§,§ 










1.170 


0.0969 


0.427 




0.715 


0.1044 


0.227 




0.576 


0.1009 


0.153 




0.412 


0.1020 


0.070 


2n-n @ 


1.316 


0.0990 


0.476 




0.997 


0.0983 


0.353 




0.831 


0.0987 


0.278 




0.730 


0.0996 


0.229 




0.603 


0.1001 


0.167 




0.560 


0.0997 


0.146 




0.425 


0.0986 


0.082 



"-LJ @ 
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0.418 0.1073 0.073 



Colloid y] 


0.712 
0.562 


0.0970 
0.1023 


0.225 
0.144 


Asymm. hard-spheres [jl4j 










0.186 


0.0744 


0.005 




0.173 


0.0758 


0.003 




0.164 


0.0788 


0.002 



TABLE I. Values of T* , r and the range of the equivalent 



square well system R for different potentials, and for different 
ranges, calculated at the liquid-gas critical point 
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